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Abstract 

We present a simple two-dimensional dynamical system reaching a sin- 
gularity in finite time decorated by accelerating oscillations due to the inter- 
play between nonlinear positive feedback and reversal in the inertia. This 
provides a fundamental equation for the dynamics of (1) stock market prices 
in the presence of nonlinear trend- followers and nonlinear value investors, 
(2) the world human population with a competition between a population- 
dependent growth rate and a nonlinear dependence on a finite carrying 
capacity and (3) the failure of a material subject to a time-varying stress 



with a competition between positive geometrical feedback on the damage 
variable and nonlinear healing. The rich fractal scaling properties of the 
dynamics are traced back to the self-similar spiral structure in phase space 
unfolding around an unstable spiral point at the origin. 

Singularities play an important role in the physics of phase transitions as well as in 
signatures of positive feedbacks in dynamical systems, with examples in the Euler equations 
of inviscid fluids in vortex collapse of systems of point vortices, in the equations of 
General Relativity coupled to a mass fleld leading to the formation of black holes 0, in 
models of micro-organisms aggregating to form fruiting bodies 0, in models of material 
failure 0, of earthquakes and of stock market crashes p. 

The continuous scale invariance usually associated with a singularity can be partially 
broken into a weaker symmetry, called discrete scale invariance (DSI), according to which 
the self-similarity holds only for integer powers of a speciflc factor A 0. The hallmark of 
this DSI is the transformation of the power law into an oscillatory singularity, with log- 
periodic oscillations decorating the overall power law acceleration towards the singularity. 
Such log-periodic power laws have been documented for many systems such as with a built- 
in geometrical hierarchy, as the result of a cascade of ultra-violet instabilities in growth 
processes and rupture, in deterministic dynamical systems, in response functions of spin 
systems with quenched disorder, in spinodal decomposition of binary mixtures in uniform 
shear flow, etc. (see and references therein). 

Here, we introduce a general dynamical mechanism for a flnite-time singularity with self- 
similar oscillatory behavior, based on the interplay between nonlinear positive feedback and 
reversal in the inertia: 

dyi 

^ = ^- 

^ = y2\y2r-' - lyilvir' . (2) 

This model is motivated by and derived from the dynamics of stock market prices, of the 



world human population and of material failure as we shall see below (see for detailed 
derivations). Our analysis of (|I], |^) offers a fundamental understanding of the observed in- 
terplay between accelerating growth and accelerating (log-periodic) oscillations previously 
documented in speculative bubbles preceding large crashes ||6|.H, the world human popula- 



tion 1 17, 16 1, and time-to-failure analysis of material rupture 0] (and references therein). 

Stock market price dynamics: the heterogeneous behavior of agents has recently been 
shown to be a crucial ingredient to account for the complexity of financial time series (see 
||T0| and references therein). Typically, "value investors" track the fundamental price of a 
given stock placing investment orders of (algebraic) size fivaiuo(i) while "technical analysts" 
use trend following strategies to place investment orders of size Qtech{t)- The balance between 
supply and demand determines the price variation from p{t) to p{t+6t) over the time interval 
6t according to In p{t+6t)—\np(t) = j^['^vaiuc(t)+Qtcch(t)] |Tl| where L > is a market depth. 
We postulate the nonlinear dependence fivaiue(^) = ~c lia[p(t)/pf] \\n[p(t)/pf]\^~^, where 
n > 1 and c > is a constant. The case n = 1 retrieves an ingredient of previous models 
|TT|,[T2|. According to textbook economics, p/ is determined by the discounted expected 
future dividends whose estimation is very sensitive to the forecast of their growth rate and 
of the interest rate, leading to large uncertainties in pf. As a consequence, a trader trying to 
track fundamental value has no incentive to react when she feels that the deviation is small 
since this deviation is more or less within the noise. Only when the departure of price from 
fundamental value becomes relatively large will the trader act. The exponent n > 1 precisely 
accounts for this effect. The second class of investors follow strategies that are positively 
related to past price moves. This can be captured by the following contribution to the order 
size: fitcch(i) = ai[lnp(t) - \np{t - 5t)] + a2[lnp(t) - \np{t - 6t)]\ \np{t) - \np{t - 5t)|'^-i 
with ai > and 02 > 0. The choice m > 1 means that trend- following strategies tend to 
under-react for small price changes and over-react for large ones. Posing yi(t) = \n[p{t)/pf], 
we expand the equation of balance between supply and demand as a Taylor series in powers 
of 6t and get 
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(Px r, o,^^ ^. dx a2{5t)^ dx dx _^ c \in-i rni/x ^3 



L 



where represents a term of the order of (St)^. This equation is a generahzation of 

the model of Pandey and Stauffer |jl3|, by allowing nonlinear positive feedback m > 1 of 
the trend-following strategies. The theory becomes critical when the "mass" term vanishes, 
i.e., when ai = L. Rescaling t and yi by a and posing y2 = dyi/dt and 7 = a~'^'^^^'>c/ L[5tY 
where a = a2(5t)''""^/iv, we obtain 

Population dynamics: the logistic equation corrects Malthus' exponential growth 
model by assuming that the population p{t) cannot growth beyond the earth carrying ca- 
pacity K: ^ = (Jop{t) [K — p{t)] where (Tq controls the amplitude of the nonlinear saturation 



term (see [jT5| and references therein). However, it is now understood that K is not a con- 
stant but increases with p{t) due to technological progress such as the use of tools and 
fire, the development of agriculture, the use of fossil fuels, fertilizers etc. as well an expan- 
sion into new habitats and the removal of limiting factors by the development of vaccines, 
pesticides, antibiotics. If K{t) grows faster than p(t), then p(t) explodes to infinity after 
a finite time creating a singularity due to the corresponding growth of the growth rate 
cr = d\np/dt, leading to ^ oc cr^ |jl6|. We now generalize it by assuming a nonlinear satu- 
ration: ^ = Q;cr|o"|''"~^ — 'y\n{p/ Koo) \ ln(p/i^'oo) 1*^"^ where a and 7 > measure the effect 
of feedback and reversal. Apart from the absolute value, the first term in the r.h.s. is the 
previous nonlinear growth of the growth rate. The novel second term favors a restoration 
of the population p{t) to the asymptotic carrying capacity K^o- The effective cumulative 
growth rate ln{p/Koo) is the natural variable to describe the attraction to Koo)- The nonlin- 
ear restoring exponent n > 1 captures the many nonlinear (often quasi-threshold) feedback 
mechanisms acting on population dynamics. Defining variables (2/1,2/2) = (a ln(p/-ft'oo), o") 
and rescaling t by a lead to ([l|J^). 

Rupture of materials with competing damage and healing: a standard model of 



damage rupture [18| consists in a rod subjected to uniaxial tension by a constant apphed 



axial force P. The undamaged cross section A{t) of the rod is assumed to be a function 



of time but is independent of the axial coordinate. The considered viscous deformation is 
assumed to be isochoric: AqLq = A{t)Lit) = constant at all times, where L{t) is the length 
of the rod. The creep strain rate ^ = = — ^^is assumed to follow Norton's law: 
^ = Co"'' where a = P/A is the rod cross section with C > and /i > Eliminating 
dec/dt leads to A^~^dA/dt = —CP^, and hence ^ = Ca^^^. Physically, this results from a 
geometrical feedback of the undamaged area on the stress and vice- versa. 

We generalize this model by adding healing as well as a strain-dependent loading: 
^ = aalcrl"^^^ — 7ec|ec|"~^. The first term in the right-hand-side is identical to previ- 
ous geometrical positive feedback for m = n — 1. We relax this correspondence and allow 
m > 1 to be arbitrary. The addition of the second term introduces the physical ingredient 
that damage can be reversible. Large deformations can enhance healing which increases 
the undamage area and thus decrease the effective stress within the material. The model is 
completed by using again Norton's law but with exponent m': ^ = Ca'"^' . Incorporating 
the constant C in a redefinition of time Ct t (with suitable redefinitions of the coefficients 
a/C —>■ a and 7/C 7), taking m' = 1 and posing (?/i,?/2) = (ccO"), we retrieve 

Analysis of the dynamical system (0,|2D for n > 1 and m > 1 with 7 > 0: When 
only the reversal term is present (i.e., the term y2\y2\"^~^ is absent), (|l|,|^) describe a non- 
linear oscillator with conserved Hamiltonian H{'y; n, 7) = ^^(^1)^ + ^l/i- trajectory 
is periodic along a constant H with period 

T{H; n, 7) = C{n) 7^ , (4) 

where C{n) is a positive number that can be explicitly calculated 0. 

When only the positive feedback is present (i.e., 7 = 0), (|l],|D leads to a finite-time 
singularity. The solution y(t;yo,^o) with initial condition yo = (1/10,^20) at time to can be 
explicitly written as: 

. m-2 

\7Tl — m — 2 m — 2 

yi{t;yo,to) -yio = sign[|/2o] — 2-m — t'^^'^'^^^o) -t)~ - (tc(z/2o) -to)~] for m 7^ 2; 

I tc — to \ ^ 



sign[|/2o] {m - 1) 



m— 1 



{tc{y2o) - 1) 



m — l 



(5) 



where the critical time interval tc(y2o)— = ^;^|2/2or~™' depends only on ?/20- Ast ^c(y2o), 
?/2 becomes either +oo or —oo depending on the sign of ?/20 for any m > 1. For 1 < m < 2, 
it is easy to show that \yi\ also becomes infinity as t — > tc(2/2o)- In contrast, for m > 2, yi 
remains finite. We thus think that m > 2 is the relevant physical regime for the financial, 
population and rupture problems discussed above. From here on, our analysis focuses on 
n > 1 and m > 2. 

Putting the nonlinear oscillation and positive feedback terms together, the dynamics of 
(0,^ is characterized in figure |l[ The two special intertwined trajectories along (solid line) 
and b~ (dashed line) connect the origin (0,0) to (+oo, +oo) and (— oo,— oo), respectively, 
and hence divide the phase space y = (yi,?/2) into two distinct basins and . The 
basin (resp. B^) corresponds to a finite-time singularity yc(yo) with ?/2c(yo) = +oo 
(resp. ?/2c(yo) = —oo) but finite yidyo) at the critical time tc(yo)- 

Starting from yo in B~^ close to the origin at to, a trajectory y(t; yo, to) spirals out with 
clockwise rotation and we count a turn each time it crosses the ?/i-axis, i.e., yi changes its 
direction of motion {dyi/dt = 0). Deep in the spiral structure, y(t;yo,to) follows approxi- 
mately the orbit of constant Hamiltonian H defined by the nonlinear oscillator but fails to 
close on itself because H is no longer conserved due to the positive feedback: 



This slowly growing nonlinear oscillator defines the first dynamical regime. 

Any trajectory starting in the first dynamical regime must eventually cross-over to the 
second one associated with a route to the singularity without any further oscillation. The 
second dynamical regime occurs ?/2 diverges (and therefore dy2/dt also diverges) while yi 
remains finite on the approach to tc(yo)- As a consequence, the reversal term 7?/i|?/i|"~^ 
in dy2/dt (0) becomes negligible close to tc(yo) and (^ is the asymptotic solution of (0,0). 
Figure ^ shows a typical time series of a y(t; yo, to) starting from yo near the origin. 



d_ 
di 



H{y;n,m,'j) = \y2 



m+l > Q _ 



(6) 
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In the basin B~^, the transition from the first to second dynamical regime occurs at the 
exit segment Ae'-^^''^-' on the positive ?/i-axis (figure |^) whose right and left end-points 
according to the forward direction of the flow are defined by the boundaries and b~ , 
respectively. In forward time, Ae'-^^'^^-' fans out rapidly over yi E {—oo, oo) as it reaches a 
singularity with y2 +00 in the second dynamical regime. Note that \yi\ can reach 00 if 
and only if yo is at an end point of Ae*^+^'^°\ i.e., on either 6+ or . Similar results hold 
for the exit segment Ae'-"^'"'"'') in B~ but with 1/2 —00. In backward time, Ae^^^'~^^ and 
/\g{-i,+o) g^j];! lYito the origin while making countable infinite many turns. These backward 
swirls of the exit segments completely define the first dynamical regime. 

The first dynamical regime exhibits remarkable scaling properties that we quantify by the 
dependence on the initial condition yo = {yio, 0) on the ?/i-axis of the following quantities: 
the number {Ntum) of turns before reaching the singularity, the exit time {te) to reach the 
exit segment into the second dynamical regime, the time interval (Atg) and the increment 
(A|/i) in the amplitude of yi over one turn. Figure |^ shows log-log plots of the scaling 
properties measured at the fc-th backward intersection of 6+ with the yi-axis starting from 
A; = at the out-most intersection (the left end-point of Ae*-^^'"''^'' in figure |l|). In order to get 
accurate and reliable estimations of these dependences and of the exponents defined below, 
we have integrated the dynamical equations using a fifth-order Runge-Kutta integration 
scheme with adjustable time step. 

The log-log dependences shown in figure |^ qualify power laws defined by 

A^turn ~ , where a > , (7) 

te ~ Vio , where 6 > , (8) 

A|/io ~ yto , where c> , (9) 

Ate ~ Vw , where d>0 . (10) 

d 

Eliminating yio between (|^) and (p!OD gives Ate ~ te ■ Since Ate is nothing but the 
difference Ate = ^e(^turn + 1) ~ '^e(^turn), this gives the discrete difference equation 
te(jVturn+AJVt„,„)-te(jVtu.-n) ^ ^| ^^ieTe AiVturn = Ntum + 1 " A^'turn = 1. This provides a discrete 

* turn 
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difference representation of the derivative dt^/dNi^T.^ which can be integrated formally to 
give A^turn ~ * , which is valid for d < b. Comparing with the relation between A^tum 

a 

and te obtained by eliminating yiQ between and i.e., A^turn ~ , we get the scaling 
relation 

a = b-d. (11) 

Since a > 0, the condition d < b is automatically verified. 

Similarly, Ayw = yio(A^turn + 1) - Z/io(A^turn) oc yl^ according to definition (|). This gives 
the differential equation dyio/dNtum ~ y^o, whose solution is A'turn ~ ^/Vio^, valid for c > 1. 
Comparing with the definition (|^, we get the second scaling relation 

a = c-l. (12) 

Since a > 0, the condition c > 1 is automatically satisfied. 

Deep in the spiral structure depicted in figure |I], in the presence of the positive feedback 
term, one rotation around the origin is not exactly closed but the failure to close, which is 
very small especially near the origin, is quantified by AyiQ shown in figure 2 which follows (^. 
We approximate the time Ate needed to make one full (almost closed) rotation by the period 
T{H) without the positive feedback term. This is essentially an adiabatic approximation in 
which the Hamiltonian H and the period T{H) are assumed to vary sufficiently slowly so 
that the local period of rotation follows adiabatically the variation of the Hamiltonian H. 
Putting together (|) and the fact that the amplitude of j/io is proportional to H~ , we get 
Ate ~ |z/io|^~, which, by comparison with (|IU]), gives 

d=!i^. (13) 

The last equation is provided by dT/dt, expressed as dT/dt = (dT/dH) x (dH/dt), where 
dT/dH is obtained by differentiating (^ and dH/dt is given by (^. Estimating dT/dH from 
(^) is consistent with the above approximation in which a full rotation along the spiral takes 
the same time as the corresponding closed orbit in absence of the positive feedback term. 



Expressing dH/dt using involves another approximation, which is similar in spirit to a 
mean-field approximation corresponding to average out the effect of the positive feedback 
term over one full rotation. In so doing, we average out the subtle positive and nega- 
tive interferences between the reversion and positive feedback terms. We replace dT/dt by 
dh.te/dte and obtain ~ H''^^^ |i/2|™^^ ~ Ate""' \yio\~ ^ , where the dependence 
in Ate is derived by replacing H by its dependence as a function of T (by inverting T{H)) 
and by identifying T and Ate- We have also used \y2\ ~ \yio\^~ ■ Taking the derivative of 



d 



Ate ~ with respect to te provides another estimation of and replacing Ate in the 



above equation by its dependence as a function of yio as defined by (|TOl) gives finally: 



a = l(n + l)(m + l)-^(3n + l) . (14) 

We find a perfect agreement between the theoretical predictions (llll), (HD) (HI)) (HD for the 
exponents a, b, c, d defined by (0)-(|lO]) with an estimation obtained from the direct numerical 
integration of the dynamical equations 0. We have also verified the independence of the 
exponents a, b, c, d with respect to the amplitude 7 of the reversal term. 

In the oscillatory regime, the growth of the amplitude of yi{t) follows a power law similar 
to (^. We obtain it by combining some of the previous scaling laws ([7|-p^. Indeed, taking 
the ratio of (|^) and (|10|) yields Ayi/Atg ~ vt^'^- Since Ayi corresponds to the growth 
of the local amplitude Ay-^ of the oscillations of yi(t) due to the positive feedback term 
over one turn of the spiral in phase space, this turn lasting Ate, we identify this scaling 
law with the equation for the growth rate of the amplitude of yi in this oscillatory 
regime: ~ ^yf^ whose solution is Ay-^{t) = ^^,^^y/b , where B is another amplitude. 

This prediction is verified accurately from our direct numerical integration of the equations 
of motion 0. We have used the scaling relations (jl^) and ([T2D leading to c + d — 1 = b. The 
time t* is a constant of integration such that B/{t*)^ = Ay-^^to), which can be interpreted 
as an apparent or "ghost" critical time, t* has no reason to be equal to te, in particular 
since the extrapolation of Ay-^{t) = (^^,^^y/b too close to t* would predict a divergence of 
yi{t). The dynamical origin of the difference between t* and te comes from the fact that t* 
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is determined by the oscillatory regime while tc is the sum of two contributions, one from 
the oscillatory regime and the other from the singular regime. 

Combining (^) with this solution for Ay-^^ {t) gives the time dependence of the local period 
Ate of the oscillation in the oscillatory regime as Atg = tfc+i —tk ~ {t* — tfc)^, where k is the 
turn index previously used. This result generahzes the log-periodic oscillation associated 
with discrete scale invariance (DSI) characterized by t^+i — ~ 1/X'' (where A > 1 is a 
preferred scaling ratio of DSI), which is recovered in the limit d/b 1~ (corresponding 
to {n oo,m 2). The dynamical system (|l],|2|) provides a mechanism for generalized 
log-periodic oscillations with, in addition, a finite number of them due to the cross-over to 
the non-oscillatory regime. We shall report elsewhere on tests of this theory on financial 
and rupture data. 
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FIGURES 




FIG. 1. Geometry of the boundaries b'^ and as well as the basins and B'^ for 
(n, m) = (3, 2.5) and 7 = 10 in phase space: the exit segments Ae^"^^'"'^^ G B^ and Ae''"^'^'^) € B~ 
are thick solid and dashed lines respectively on the 7/1-axis; the squares and diamonds are the initial 
condition and exit point of the trajectory in Figure |^. 
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FIG. 2. A trajectory in B+ starting from initial condition (-2.32878 x 10^^,1.71083 x 10"^) 
at t = —500, going through the exit point (0.5, 0) at t = 0, to the critical point (0.768899, +oo) at 
tc = 2.67284: Solid and dashed lines correspond to yi and y2, respectively; square, diamond and 
triangle are the initial condition, the point in the exit segment Ae^^^'^*^), and the critical point 
along the trajectory with open and filled symbols corresponding to yi and ?/2- 
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FIG. 3. Scaling laws associated with the self-similar properties of the nonlinear oscillatory 
regime as a function of initial conditions at turn points for (n, m) = (3, 2.5) and 7 = 10. 
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